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VIII. On the Influence of Geological Changes on the Earth's Axis of Botation. By 
Geoege H. Darwin, M.A., Felloio of Trinity College^ Cambridge. Communicated 
by Professor J. C. Adams. 

Eeeeived October 13, — Read November 23, 1876. 

The subject of the fixity or mobility of the earth's axis of rotation in that body, and 
the possibility of variations in the obliquity of the ecliptic, have from time to time 
attracted the notice of mathematicians and geologists. The latter look anxiously for 
some grand cause capable of producing such an enormous effect as the glacial period. 
Impressed by the magnitude of the phenomenon, several geologists have postulated a 
change of many degrees in the obliquity of the ecliptic and a wide variability in the 
position of the poles on the earth ; and this, again, they have sought to refer back to 
the upheaval and subsidence of continents. 

Mr. John Evans, F.R.S., the late President of the Geological Society, in an address 
delivered to that Society, has recurred to this subject at considerable length. After 
describing a system of geological upheaval and subsidence, evidently designed to produce 
a maximum effect in shifting the polar axis, he asks : — " Would not such a modification 
of form bring the axis of figure about 15° or 20° south of the present, and on the 
meridian of Greenwich — that is to say, midway between Greenland and Spitzbergen % 
and would not, eventually, the axis of rotation correspond in position with the axis of 
figure % 

" If the answer to these questions is in the affirmative, then I think it must be conceded 
that even minor elevations within the tropics would produce effects corresponding to 
their magnitude, and also that it is unsafe to assume that the geographical position of 
the poles has been persistent throughout all geological time. 5 '* 

On the few occasions on which this subject has been referred to by mathematicians, 
the adequacy of geological changes to produce effects of such amount has been denied. 
Amongst others, the Astronomer Royal and Sir William Thomson have written briefly 
on the subjectf , but, as far as I know, the subject has not hitherto been treated at much 
length. 

The following paper is an attempt to answer the questions raised by Mr. Evans ; but 
as I have devoted a section to the determination of the form of continent and sea which 
would produce a maximum effect in shifting the polar axis, I have not taken into consi- 
deration the configuration proposed by him. 

The general plan of this paper is to discuss the following problems : — 

* Quart. Journ. Geol. Soc. 1876, xxxii, Proc. p. 108. f In papers referred to below. 
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First. The precession and nutation of a body slowly changing its shape from internal 
causes, with especial reference to secular alterations in the obliquity of the ecliptic. 

Second. The changes in the position of th e earth's axis of symmetry, caused by any 
deformations of small amount. 

Third. The modifications introduced by various suppositions as to the nature of the 
internal changes accompanying the deformations. 

In making numerical application of the results of the previous discussions to the case 
of the earth, it has of course been necessary to betake one's self to geological evidence ; 
but the vagueness of that evidence has precluded any great precision in the results. 

In conclusion I must mention that, since this paper has been in manuscript, Sir 
William Thomson, in his Address to the Mathematical Section of the British Association 
at Glasgow, has expressed his opinion on this same subject. He there shortly states 
results in the main identical with mine, but without indicating how they were arrived at» 

The great interest which this subject has recently been exciting both in England and 
America, coupled with the fact that several of my results are not referred to by Sir 
William Thomson, induces me to persist in offering my work to the Eoyal Society. 

I. PEECESSION OP A SPHEEOID SLOWLY CHANGING ITS SHAPE. 

I begin the investigation by discussing the precession and nutation of an ellipsoid of 
revolution slowly and uniformly changing its shape. The changes are only supposed 
to continue for such a time, that the total changes in the principal moments of inertia 
are small compared to the difference between the greatest and least moments of inertia 
of the ellipsoid in its initial state. 

For brevity, I speak of the ellipsoid as the earth ; and shall omit some parts of the 
investigation, which are irrelevant to the problem under discussion. 

The changes are supposed to proceed from internal causes, and to be any whatever ; 
and in the application made they will be supposed to go on with a uniform velocity. 

1. The Equations of Motion. 

M. Liouville has given the equations of motion about a point of a body which is 
slowly changing its shape from internal causes*; these equations, he says, are only 
applicable to the case of the point being fixed or moving uniformly in a straight line. 
They may, however, be extended to the motion of the earth about its centre of inertia, 
because the centrifugal force due to the orbital motion and the unequal orbital motion 
will not add any thing to the moments of the impressed forces. These equations are, 
in fact, an extension of Euler's equations for the motion of a rigid body, which are 
ordinarily applied to the precessional problem. To make them intelligible I reproduce 
the following from Mr. Eouth's ' Eigid Dynamics ' f, where the proof is given more 
succinctly than in the original : — 

* Liotjville's Journ. Math. 2 me serie, t. iii. 1858, p. 1. 

f Page 150, edit, of 1860, but omitted in later editions. L, M } N" are the couples of the impressed forces 
about the axes. 
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" Let x\ y> z be the coordinates of any particle of mass m at the time £, referred to 
axes fixed in space. Then we have the equation of motion 

and two similar equations, 

"Let 7 ^ / dy doc\ 

ll *=^ m \ x Tt-yjt) . • . . (2) 

with similar expressions for A 1? h 2 . 
Then the equation (1) becomes 

" Let the motion be referred to three rectangular axes Ox\ Oy\ Oz 1 moving in any 
manner about the origin O. Let a, /3, y be the angles these three axes make with the 
fixed axis of z. Now h 3 is the sum of the products of the mass of each particle into 
twice the projection on the plane of xy of the area of the surface traced out by the radius 
vector of that particle drawn from the origin. Let h\, h' n A', be the corresponding 
< areas ' described on the planes y ! z ! , z f x f , x f y f respectively. Then by a known theorem 
proved in Geometry of Three Dimensions, the sum of the projections of A"/, A' 2 , A' 3 on xy 
is equal to h 3 ; 

.-. A 8 =A' 1 cos ct+h'i cos 0+A' 3 cos y ....... (4) 

" Since the fixed axes are quite arbitrary, let them be taken so that the moving axes 
are passing through them at the time t. Then 

fl 1 =zfi l ^ fi 2 =zfi 2 , .fi 3 zzzfi 3 ', 

and by the same reasoning, as in Arts. 114 and 115, we can deduce from equation 
(4) that 

tiA + MA (5) 



dh 8 ___ dh ! 3 
dt dt ' 



where 0„ 6 29 6 3 are the angular velocities of the axes with reference to themselves. 
Hence the equations of motion of the system become 



dh ^ 

-^ li 2 3 + h 3 6 2 = L, 

dh c 



J 



> • • •■ (6) 



" These equations may be put under another form which is more convenient. Let 
d, y\ %' be the coordinates of the particle m referred to the moving axes, and let 
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" Since the fixed axes coincide with these at the time t 9 we have w=tf 9 y=y\ and by 

Jl\.x t/» J L\jt, 

d l- d l+b x b z 



and by similar reasoning 



dt dt _j 

.-. //^IF. + C^-E^-D^*; 



" Hence the general equation of motion becomes 

j t (Cfl,-E^-Dfl a +H,)+F(fll-fl?)+(B-A)flA+EflA-DflA+a,H 2 -fl 8 H I =N . (7) 

and two similar equations. 

"Let the moving axes be so chosen as to coincide with the principal axes at the 
time t. Then D:=0, E=0, F=0, and the equations become," 

and two similar equations; where a 1? X 29 x 3 (replacing the A, B, C of Mr. Routh) are 
the three principal moments of inertia, and are functions of the time. 

In order to apply these equations to the present problem, we must consider the 
meaning of the quantities Q i9 2 , 3 . A system of particles may be made to pass from 
any one configuration to any other by means of the rotation of the system as a whole 
about any axis through any angle, and a subsequent displacement of every particle in a 
straight line to its ultimate position. Of all the axes and all the angles about and 
through which the preliminary rotation may be made, there is one such that the sum 
of the squares of the subsequent paths is a minimum. By analogy with the method of 
least squares this rotation may be said to be that which most nearly represents the 
passage of the system from one configuration to the other. If the two configurations 
differ by little from one another, and if the best representative rotation be such that 
the curvilinear path of any particle is large compared to its subsequent straight path, 
the system may be said to be rotating as a rigid body, and at the same time slowly 
changing its shape. Now this is the case we have to consider in a slow distortion of 
the earth. 

Divide the time into a number of equal small intervals r, and in the first interval let 
the earth be rigid, and let each pair of its principal axes rotate about the third (with 
angular velocities double those with which they actually rotate). At the end of that 
interval suppose that each pair has rotated about the third through angles 2^r, 2a> 2 r 9 

* A, B, C, D, E, F are, as usual, the moments and products of inertia. 
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2co 3 r. Then reduce the earth to rest, and during the next interval let the matter 
constituting the earth flow (with velocities double those with which it actually flows) so 
that the pairs of principal axes have, at the end of the interval, rotated with respect to 
the third ones through the angles — 2ar, — 2/3r, — 2yr. Lastly let 25 x r, 23 2 r, 2$ 3 r be 
the rotations of each pair of axes about the third by which they could have been 
brought directly from their initial to their final positions in the time 2r. 
Therefore, by the principle of superposition of small motions, 

Now supposing these two processes to go on simultaneously with their actual velocities, 
instead of in alternate intervals of time with double velocities, it is clear that 6 U 2 , $ 3 
are " the angular velocities of the axes with reference to themselves"; a Xi a/ 2 , co 3 are the 
component angular velocities of the earth considered as a rigid body; and —a, -—0, — y 
are the component angular velocities of the principal axes relatively to the earth, arising 
from the supposing continuous distortion of that body. 

With respect to the other quantities involved in the equations of motion: — 
Let C, A be the principal moments of inertia of the earth initially when t is zero ; 
and at any time t, let 

\=A + a£, \ 2 =A+bt, X 3 =C+c£. 

We here suppose that the changes in the earth are so slow that terms depending on 
higher powers of t may be neglected. 

Lastly the quantities H n H 2 , H 3 are respectively twice the areas conserved on the 
planes of #A, 6 A, && by the motion of the earth relative to these axes. If the earth 
were rigid, they would all be zero, because there would be no motion relative to the 
principal axes : thus <y M a) 2 , a> 3 do not enter into these quantities. Now the motion which 
does take place may be analyzed into two parts. Divide the time into a number of 
equal small elements r, and in the first of them let the matter constituting the earth 
flow (with a velocity double that with which it actually flows) ; this motion will 
conserve double-areas on the planes of 2 3 , QJ& 19 &$„ which we may call 2^ 1 r,2^ 2 r, 2$J 3 r. 
In the next interval of time let each pair of axes rotate round the third (with angular 
velocities double those with which they actually rotate), so that at the end of the 
interval they have turned through the angles — 2ar, — -20r, — 2yr. Now since during 
this second interval the axes have rotated in a negative direction through the solid, 
therefore the solid has rotated in a positive direction with reference to the axes. 
Remembering then that X„ A 2 , X 3 are the principal moments of inertia, the double-areas 
conserved on the three planes in this second interval are 2x x ar, 2X 2 j3r, 2x 3 yr. Hence if 
2H 1 r, 2H 2 r, 2H 3 r be the double areas conserved in this double interval of time, we have 
2H 1 r=2l I r+2?^r, 2H 2 r==21 2 r+2x 2 £r, 2H 3 r=z21 3 r+2x 3 yr. 

Therefore if we now suppose the two processes to go on simultaneously with their actual 
velocities, instead of in alternate elements of time with double velocities, and if we 
substitute for x„ X 2 , X 3 their values in terms of A, a, t f &c, we get 
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H 1 =(A+a«)«+^ 1 , H 2 =(A+W)0+1 2 , H,=(C+c*)y+= 

where Jh, ^ 2 , H 3 , denote those parts of the double areas conserved, which depend only 
on the internal motions accompanying the change of shape. 

Then, if the changes proceed with uniform velocity, a, a, ^. &c. are all constant. 

Corresponding also to the equations of motion are the geometrical equations 



dt 



^aCOsp-fflxSiiKp ; 



d^f , A A A 

-n sm9= — t^ cos <p+9 2 sin <p ; 

d<p . dfy A A 
dt+dt C0B(> = ^ 




In figure 1, A, B, C are the axes, about which the moments of inertia are X l5 X 2 , ^ 3 ; 
XY is the ecliptic ; and the meaning of the other symbols is sufficiently indicated. 

Substituting, then, for the various symbols in the original equations of motion, it will 
be found that 



A-^ — (A' — G)6>2*> 3 =li— £<a-^r— -(b — c)a> 2 a) B +byt 
and two similar equations*. 



ej 2 ^C^ej B >--^ l —{AJy+%}cJ 2 



ym, 



d&Y 



Now the terms on the right-hand side are always very small compared to K~n , because 

the time will not run on until they have become large ; hence approximate values may 
be substituted therein. 



* The A is written A' in two places, where it may he taken to stand for B ; and then the other equations 
may be found by cyclic changes of letters and suffixes. 
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Let the angular velocity of rotation of the earth be — n, and let II cosec Q be the pre- 
cession of the equinoxes ; then in the small terms the following substitutions may be 

made :— 

<p = — (nt+e), #!= — Hcos(nt-\-e), eo 2 =— XIsin(^-f-£) 5 ^3=— n, 

and the e may be omitted for brevity. 

Further, N (depending on the attractions of the sun and moon) is very small ; and a 
consideration of the third equation of motion shows that, when integrated, it leads to 
^3=—^+ terms, which are very small during the limited period under consideration. 
And if these terms were substituted on the left-hand side of the two former equations, 

C— A 
they would be still further diluted by multiplication by the small quantities — -^— 

and II. Hence the third equation may be neglected, and in the two former equations 
—n may be substituted for & s . 

C— A 

Making these substitutions, then, and writing ^ for — r~ ^» the equations become 



doo l 
If 



dcog 
~di 



L Ilnf . , by) , . ■ , Ila , irr ( , p?o) . , 

^ 2 ==-r —j- <a — b+c—~^>tsmnt+-j- cos nt+Yliy+j^ysmnt 

-£{ A - 3 }-X^ + y) ; 

M TIn( , . . ay) , , , lib . , —( . P? 3 \ , 

+^ 1 ==^~.~^^a+b+c~^~>£^ 



+ a \ A 



Then neglecting ~ compared to unity, putting C + ct=A in the small terms, and only 
retaining the more important terms, 

da), L Tin r , , , , . , , Ila # 1 T-r f i S?s ) • j. o 

-^ — ^ 2 =^— j{a- b + c}t sm nt+j^ cob nt+U<y+ ^>Bmnt—n(3 9 



^+p®i~A+~x\— &+b+c}t cos nt+--£ amnt—IlW+^lcoant+m*. 



3{„„„' M 4- f /w^* 



These are the required equations of motion, and in integrating them they may be 
treated as linear, 

* If we wish to treat a, a, ffi v &c. as variable, we have only to add to the right-hand sides of these equations 

l lt cos nt ^--~ Si 3 and 5 1 sin nt d ^~l-^ respectively. If we put L=M=0 and neglect %, p| 2J %, 
A dt A dt A dt A dt * J 

these equations will "be found to be identical with the equations (2) given by Sir "W. Thomson in App. C. I 
had not noticed until it was pointed out by him, how nearly applicable my equations were to the case of 
varying velocities of distortion.-— April 26, 1877. 
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2. Inequalities independent of the Impressed Forces. 

First, then, suppose that L=M=0. 

Integrate the equations, and neglect ^ compared with w, and we have 

TT XT TT / wc* \ 

+-( a+~r ) + r cos for —Cj smft»r, 

TT I I' '"i " ' I* / 5ft* \ 

® % — A (— a + b +c) £ sin nt+j^ (c— a) cos rc#— —fy+^j sin w£ 

+ -( /3+-V" 2 j — F sin^— Gcos ^1 

The last two terms in co l and oj 2 represent the complementary function ; and the values 
of F and G must be determined from the initial conditions. 

Now it will be shown later that a, 0, y are comparable with tt^t^ n ^t ? p —~A ; hence 

the terms in the second lines are much more important than those in the first. Thus 
in determining the values of F and G we may neglect the first lines. 

Now, initially, the instantaneous axis coincides with the axis of greatest moment of 
inertia; so that when t=0, ^=^=0, and therefore 



f =-^+i). g =#+ 



£jj? 



The terms in F and G represent an inequality of 306 days period. 

3. The Inequality of 306-Days Period*. 

I have worked out the value of $? x and % 2 in two supposed cases of elevation, under 
certain suppositions as to the nature of the internal movements of the earth. In one of 

* I have thought it necessary to discuss this inequality fully, both on account of its intrinsic interest, and 
because it has been referred to by the Astronomer Royal and Sir "William Thomson. 

The former says Q Athenaeum/ Sept. 22, 1860) :— a Fow, let us suppose the earth not absolutely rigid, but 
that there is susceptibility to change of form, either from that degree of yielding or fracture to which most 
solid substances are liable, or from the hydrostatic pressure of internal fluid. This, as I conceive, puts an end 
to all supposition of change of axis. The first day's whirl would again make the axis of rotation to be a prin- 
cipal axis, and the position of the axis would then be permanent." 

Eut Sir Geokge Airy is here speaking of the effect of the elevation of a mountain mass in about latitude 45°, 
by something like a gaseous explosion. This supposition is not at all in accordance with the belief of geologists, 
whereas a gradual elevation is so. 

Sir W. Tiiomsox, on the other hand, says (Trans. GeoL Soc. Glasgow, 1874, vol. xiv. p. 312) :— " In the 
present condition of the earth, any change in the axis of rotation could not be permanent, because the instan- 
taneous axis would travel round the principal axis of the solid in a period of 296 days. .... In very early 
geologic ages, if we suppose the earth to have been plastic, the yielding of the surface might have made the new 
axis a principal axis. Eut certain it is that the earth at present is so rigid that no such change is possible/' 
And he adds that practical rigidity has prevailed throughout geologic history. 
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them I found fx— fTp and H x = ; and in the other ^ =^ 2 = 0. In order not to inter- 
rupt the thread of the argument, the calculation is given in Appendix A; it will also 
be more intelligible after the latter part of this paper has been read. In the general 
case the same kind of proportion will subsist between $h and Aa, ^ 2 and Aj3, and we 
may therefore, without serious error, neglect the former compared with the latter. 
Thus, as far as concerns the present inequality, 

TlCt Tip 

#!=:-— (1-— cos^£)— - — sin ffit, 

co. 2 ~ — sin ^,£+ — (J- — cos i^v» 

On account of this inequality the greatest angular distance (in radians) of the instan- 

taneous axis from the pole is - Wa ' + p -. It will appear from the latter part of this 

paper that, if the elevation of a large continent proceeds at the rate of two feet in a 
century, s/^TW may be about T £o" P er annum, and p is 360° in 306 days ; whence it 
follows that the greatest angle made by the instantaneous axis with the axis of figure is 
comparable with 3^-/, a quantity beyond the power of observation. On the score of 
these terms the instantaneous axis will therefore remain sensibly coincident with the 

axis of figure. 

They will, moreover, produce no secular alteration on the obliquity of the ecliptic, nor 

in the precession, because they will appear as periodic in j t and j t sin 0, with arguments 

n and n^zp. 

Now although this inequality is so small, it nevertheless is of interest 

If we map, on a tangent plane to the earth at its initial pole, the relative motion of 

the instantaneous axis and the pole of figure, we get, as the equation to the curve, 

%=- (1 — cos put) — - - sin ysl, 
^=-sin^+ (1 — cos f*t). 
If t be eliminated from these equations, we get 



(?~i) + (y + 



V-l f* 2 



Thus the relative motion is a circle, passing through the origin, and touching a line 
inclined to the axis of y at an angle arc tan |. Therefore the instantaneous axis describes 
a circle passing through the pole of figure every 306th day; and this circle touches the 
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meridian, along which the axis of figure is travelling with uniform velocity, in conse- 
quence of the geological deformation of the earth. 

The motion of the instantaneous axis in the earth is a prolate cycloid. 

4. Adjustments to a Form of Equilibrium. 

If the earth were a viscous fluid there is no doubt but that the pole of figure would 
tend to displace itself towards the instantaneous axis, whose mean position would be the 
centre of the circle above referred to. 

But Sir William Thomson has shown* that the earth is sensibly rigid; and in any 
case the earth is not a viscous fluid, properly so called, although it may be slightly 
plastic. 

M. Teesca has shown that all solids are plastic under sufficiently great stresses, but 
that, until a certain magnitude of stress is reached, the solid refuses to flowf . Now in 
the case of a very small inequality like this, the stresses introduced by the want of coin- 
cidence of the instantaneous axis with the axis of figure are very small, even when at 
their maximum ; and every 306th day they are zero. It seems, therefore, extremely 
improbable that the stresses can be great enough to bring the earth into what M. Teesca 
calls the state of fluidity ; and therefore it is unlikely that there can be any adaptation 
of the earth's form to a new form of equilibrium in consequence thereof. 

In all the other inequalities introduced, whether arising from the first three terms 
above given in co i and &> 2 , or arising from the impressed forces, to be treated hereafter, 
the centre of the positions of the instantaneous axis is coincident with the pole of 
figure, and therefore there can hardly be any adaptation of figure eccentric to the axis 
of greatest moment to balance the stresses introduced by centrifugal force. 

It would appear probable that, whilst a geological change is taking place, the earth 
is practically rigid for long periods. But as the earth comes to depart more and more 
from a form of equilibrium, the stresses due to the mutual gravitation of the parts, and 
to the rotation, increase gradually, until they are sufficiently great to cause the solid 
matter to flow. A rough kind of adjustment to a form of equilibrium would then take 
place. The existence of continents, however, shows that this adjustment does not take 
place by the subsidence of the upheaved part ; and as this adaptation of form would 
be produced by an entirely different cause from that to which the upheaval was due, that 
upheaval would probably persist independently of the approximate adoption of a new 
form of equilibrium by the earth. 

M. Teesca's experiments on the punching of metals would lead one to believe that 
the change would take place somewhat suddenly, and would, in fact, be by an earth- 
quake, or a succession of earthquakes. On each of these occasions the tendency would 

* In his Address to the British Association, 1876, he states that the argument derived from precession 
(Thomson and Tait's Nat. Phil. p. 691) is fallacious ; he adduces, however, a number of cogent arguments on 
this point. 

t " Sur Tecoulement des Corps Solides," Mem. des Sav. torn, xviii. 
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be to adjust the form to one of equilibrium about the instantaneous axis. Now the 
principal axis X 3 has (in consequence of the postulated deformation) been travelling 

along the meridian in longitude tt-{- arc tan g, measured from the plane containing 

X x and X 3 . 

The earthquake will take place when, to the stresses due to mutual gravitation, are 
superadded the maximum stresses due to centrifugal force ; that is to say, when the 
instantaneous axis is at its greatest distance from x 3 , the axis of greatest moment of 
inertia. At the instant of the earthquake the principal axis will be moved towards the 
position of the instantaneous axis. And as the circle described by the instantaneous 
axis touches the meridian of displacement of the principal axis, therefore the principal 
axis will be carried by the adjustment towards the centre of the circle described by the 
instantaneous axis, and therefore perpendicular to the meridian of displacement. 

Thus, if the adjusting earthquakes take place at long intervals, the motion of the 
principal axis will not deviate sensibly from continuing along the meridian, along which 
it would travel in consequence merely of the geological deformation. If, however, the 
adjustments are frequent, the path of X 3 will diverge sensibly from the meridian along 
which it started. If the readjustments become infinitely frequent and infinitely small, 
there is a continuous flow of the matter of the earth, which is always seeking to bring 
back the earth's figure to one of equilibrium, from which figure it is also supposed to 
be continuously departing under the action of internal forces. In this state the earth 
may be considered as formed of a stiff viscous fluid. 

According to these ideas, at each adjustment x L , X 2 , X 3 will be suddenly reduced to 
nearly their primitive values, A, A, C ; but a, /3, y depend on the rate of accession and 
diminution of matter at various parts of the earth, and remain constant. The only 
effect, then, is that each adjusting earthquake must be taken as a new epoch. 

As far as I can see, it seems quite possible that the earth may be sensibly rigid to the 
tidally deforming influences of the sun and moon, and yet may bring itself back from 
any considerable departure from a form of equilibrium to approximately that form. It 
therefore seems worth while to consider the case of the adjustments being continuous, 
whilst the deformation is also continuous. 

5. Adjustments to the Form of Equilibrium continuous. 

I therefore propose to consider geometrically, but not dynamically, the paths of the 
instantaneous axis, and of the principal axis, where the earth is viscous and continuously 
deformed by internal forces. It is supposed that the velocities of flow of the matter of 
the earth are so small that inertia may be neglected, and that the displacements are so 
small that the principle of the superposition of small motions is applicable. 

As before the paths of the instantaneous and principal axes may be mapped on a 
tangent plane to the spheroid, at the extremity of the primitive pole, the mean radius 
of the spheroid being taken as unity, 

jQ K Zi 
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In consequence of the continuous deformation, the principal axis travels with a linear 
velocity (on the map) -\/a 2 +(3 2 along the meridian of longitude arc tan g. Take this 

meridian as axis of a?, and measure y, so that the angular velocity p is from x towards y, 
and call *Ja 2 +fi% u. 

Then the principal axis A 3 moves along the axis of x with a uniform linear velocity 
— % and, from dynamical principles, the instantaneous axis I moves round the instanta- 
neous position of A 3 with a uniform angular velocity p. 

But because of the earth's viscosity, X 3 always tends to approach I. The stresses 
introduced in the earth by the want of coincidence of X 3 with I vary as A 3 I. Also the 
amount of flow of a viscous fluid, in a small interval of time, varies jointly as that interval 
and the stress. Hence the linear velocity (on the map), with which X 3 approaches I, 
varies as X 3 I (equal to r suppose). Let this velocity be w, where v depends on the vis- 
cosity of the earth, diminishing as the viscosity increases. 

Thus the principal axis describes a sort of curve of pursuit on the map ; it is animated 
with a constant velocity —u parallel to x 9 and with a velocity vr towards I, which rotates 
round it with a uniform angular velocity ^. 

The motion of I, relative to X 3 , is that of a point moving with a constant velocity u 
parallel to #, rotating round a fixed point with a constant angular velocity p, and moving 
towards that point with a velocity vr. 

Let £, 7} be the relative coordinates of I with respect to X 3 , and #, y the coordinates 
of X 3 . Then the differential equations which give the above motions are : — 

^ — u—v% — pr t , (1) 

J=-^+^> • (2) 

^=- w + i£, (3) 

$=« w 

If (1) and (2) be integrated, and the constants determined so that, when £=0, g=yj=0 
(which expresses that initially X 3 and I are coincident), it will be found that 



u 
v 2 + ft 2 



< v(l — e~ vt cos (jbt) -\-pe~~ vt sin pi >, 

^=^t — 2")K1 —e~ vt cos {hi)—ve~ vt sin ysb\. 

These give the path of I relative to X 3 . It may be seen to be a spiral curve dimi- 
nishing with more or less rapidity, according as the earth is less or more viscous. If 
y = 0, it becomes the circle, found above from the dynamical equations. 
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Substitute in (3) and (4) for «£ and tj; integrate, and determine the constants, so that 
r hen ^=0, x=y=0. It will then be found that 



u 



y 

p(v 2 —(jtj 2 ) -+ p*t-\- ttt% \ ~~ (v* — p*)^"* cos (d-\-Zpve~ vi sin pb\ 



v^ + j"' L v r / * r • v^-j-p 
u 



J V 2 •+- jU, 2 



2p 2 +p#+ v 2TT2{2p^~^ cos (jut + (f 2 — p?)e" vt sin [d 



These give the path of X 3 on the map. It may be seen to be a cycloidal curve, in 
which the radius of the rolling circle diminishes with more or less rapidity, according 
as the earth is less or more viscous. 

After some time e~ vi becomes very small, and the motion is steady ; and then 

£= wv rj=z~-- v t or I is fixed, relatively to X 3 , at a distance —?=-== from it, and on 

h ^ 2 + v 2 ' v 2 + ju. 2 ' ' J 3 ' VV + ju 2 

the meridian, measured from the axis ef * in longitude arc tan f . This point is the 

centre of the above-mentioned spiral curve. 

If v be very small (or the earth nearly rigid) this meridian differs by little from the 
axis of y. But it may be that v is so small that e" vt has not time to become insensible 
before the geological changes cease. This case corresponds very nearly to the hypo- 
thesis, in the last section, of adjusting earthquakes. 

If the earth be very mobile, or v large, £=-, ?y=0. 

Again, with respect to the path of X 3 , when the motion has become steady, 

__ wv(v 2 —- jU, 2 ) \f?U 

X— iaqp^r^— ^qp^« h 

2^V 2 jIMV 

and eliminating t, vx^^y — — uv 2 . 

That is to say, when the motion is steady, X 3 moves parallel to meridian longitude 

5T— arc tan -, and distant from it .— ^- on the negative side. This straight line is the 

degraded form of the above-mentioned cycloidal curve. 

If the earth is nearly rigid this path does not differ sensibly from the axis of x ; if 
very mobile, it is nearly perpendicular to the axis of #, and a long way from the origin. 
In this last case the solution becomes nugatory, except as showing that the very small 
inequality of 306 days would be capable of disturbing and quite altering the path of the 
principal axis, as arising merely from geological changes on the surface of the earth. 

In the case contemplated by the Astronomer Royal, where the elevation is explosive, 
u must be put equal to zero, and the constants of integration so determined, that when 
£=0, |=E suppose, and q=x=y=Q. It will then be found that when the agitation has 
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subsided, w=K, ~, ^=B,-, or the pole of figure will have taken up a position on one 

side of the meridian, along which it was initially propelled by the explosion. 

It thus seems probable that during the consolidation of the earth there was a great 
instability in the position of the principal axis, and therefore also of the axis of rota- 
tion which followed it. 

6. Secular alteration in the obliquity of the Ecliptic, resulting from terms independent 

of the Impressed Forces. 

To return to the main line of the inquiry : — If the values of co 1 and oo^ found in sec. 2, 

be substituted in the geometrical equations for ~r and -jr sin 6 (see sec. 1), a number 

of periodic terms will arise, and these terms have diurnal and semidiurnal periods, but 
their amplitudes are so small that they have no practical interest. 

The only thing which concerns us is to inquire whether there can be any secular 
change in the obliquity of the ecliptic. 

Select, then, only terms in sin nt in a/ 15 and in cosnt in a> 2 , and substitute in the geo- 

metrical equation -r.= ---a l mint -\-oo 2 cos nt, and reject periodic terms. It will then be 
found that 

fr=-^( a + b - 2c )- 

7. Terms dependent on the Impressed Forces. 

It now remains to consider the effect of the impressed forces on the precession and 
obliquity of the ecliptic. 

The equations of motion are reduced to 

dco l C— A L 

^r + ~AT^ 2 --A' 

doo 2 C — A __M 

dco 3 __N 

It ~A' 

If we write for L, M, N, L+c!L, M+&M, cJN, and indicate by L and M the couples 
caused by the attractions of the sun and moon on the protuberant parts of the earth 
before it has begun to change its shape, then L and M only cause the ordinary precession 
and nutations. For the present problem it is therefore only necessary to consider the 
effects of ciL, £M, c>N, which arise from the change of shape of the earth. 

It follows, from the same arguments that were used in sec. 1, that the change in the 
earth's angular velocity of rotation due to <5N will only have a very small effect on oo x 
and a 2 ; so that, as far as is now important, co 3 may be put equal to — -n in the first two 
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equations, which may then be written 



dco l &L 

dt 



+^,=x- 



Now SL and SM are the changes in L and M, when A + a#, A+b£, C+c£ are written 
for A, A, and C respectively. If SL, SM be thus formed, and the equations integrated, 
it will be found that the principal terms, arising from the sun's attraction, are nine 

both in ~r. and -rr sin Q ; the same number of terms arise in the precession and nutation 

with respect to the plane of the lunar orbit, and these would have to be referred to the 
ecliptic. Sixteen out of the eighteen terms represent, however, only very small nutations, 
and the only terms of any interest are those which give rise to a secular change in the 
obliquity of the ecliptic. These terms may be picked out without reproducing the long 
calculation above referred to, for they arise entirely out of the constant couple acting 
about the equinoctial line, which gives rise to the uniform precession. 

Now this constant couple is Clin; whence L=CEfesin^, M= — CEtacos^. And 

since II involves — ~ , therefore 

C 

SL = — CUn q3t t $m ^ SM = — CEto q^\ * cos n ^- 

If these be substituted in the equations of motion and the equations integrated, and 
only terms in sin nt in u i9 and those in cos nt in cj 29 be retained, we get 

n b— c . n c— a 

__._„__ __ - sin nt, co =-—-~ ^ — x cos nt. 

1 n (J — A ' " n C — A 

Substituting in the geometrical equation ^=— ^ sin nt-\-co 2 cos nt, and rejecting 
periodic terms, 

dQ II a + b — 2c 

Jt~2ii ~~C^A " 

8. General result with respect to the Obliquity of the Ecliptic. 
It was found in sec. 6 that the secular rate of change of 6, as due to the internal 

changes in the earth, was ■— — , a+ ~ c . Since C—A is small compared to A, this term 

2n A 

is small compared with the term found at the end of sec. 7. Hence, finally, taking all 
the terms together, we get the approximate result, 

<Z0__n a + b--2c 
dt~2n C-A ' 

and for small changes in the obliquity, insufficient to materially affect II, 

. Ill a-f-b— 2c, 
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This equation has been obtained on the supposition that the change in the earth's 
form never becomes so great that at, ht, ct exceed small fractions of C— A ; a condition 
which is satisfied in the case of such geological changes as those of which we have any 
cognizance at present. 

It will appear from a comparison with results given hereafter, that — >- , __ /" • t cannot 

IX 

ever exceed two or three degrees ; and since — is a very small fraction, it follows that 

the obliquity of the ecliptic must have remained sensibly constant throughout geological 
history*. Also the instantaneous axis of rotation must always have remained sensibly 
coincident with the principal axis of figure, however the latter may have wandered in the 
earth's body. 

It has hitherto been assumed that the change of form and the angular velocities of 
the principal axes in the earth's body are uniform. But the preceding investigation 
shows clearly that no material change would be brought about by supposing the 
changes to proceed with varying velocities. This being so, dynamical considerations 
may be dismissed henceforth ; and accordingly the next part of this paper will be devoted 
to the kinematical question, as to the change in position of the earth's axis of figure as 
due to geological changes. 

The various assumptions made above will incidentally be justified in the course of 
the work. 

For some remarks of Sir William Thomson on this part of the paper see Appendix C. 

* During the Glacial Period there must have been heavy ice-caps on one or both poles of the earth. The 
above equation will give the disturbance of the obliquity of the ecliptic produced thereby. 

I will take what I believe is the most extreme view held by any geologist. Mr. Belt is of opinion that an 
enormous ice-sheet, which was thickest in about lat. 70° "N. and $., descended from both poles down to lat. 45° ; 
the amount of ice was so great that the sea stood some 2000 feet lower than now throughout the unfrozen 
regions between lat. 45° N. and S. 

Suppose that the whole of this equatorial region was sea, and that the water contained in 2000 feet of depth 
of this sea was gradually piled on the polar regions in the form of ice. Then the effect in diminishing C and 
increasing A cannot be so great as if the whole of this mass were subtracted actually from the equator and 
piled actually on the poles. This latter supposition will then give a superior limit to the amount of alteration 
in the obliquity of the ecliptic. I have calculated this alteration by means of the above formula, taking the 
numerical data used later in this paper, and taking the specific gravity of water to that of surface-rock as 4 to 11. 
I find, then, that the superior limit to the increase of the obliquity of the ecliptic would be -00045"; that is 
to say, the position of the arctic circle cannot have been shifted so much as half an inch. And this is an accu- 
mulated effect, and the matter is distributed in the most favourable manner possible. 

In this case the amount of matter displaced is enormous, and is placed in the most favourable position for 
affecting the obliquity; hence, a fortiori, geological changes in the earth cannot have sensibly affected the obliquity. 

But although this equation leads to no startling results in the geological history of the earth, I hope to show 
in a future paper that it may have some bearing on the very remote history of the earth and of the other planets 
(see a paper "On a Suggested Explanation of the Obliquity of Planets to their Orbits," Phil. Mag. March 1877). — - 
[In consequence of a mistake in the work it was erroneously stated in the abstract of this paper in the ' Pro- 
ceedings ? that the change in the position of the arctic circles might amount to 3 inches, instead of to half an 
inch. — Added August 18, 1877.] 
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II. THE PKIjSFCIPAL AXES OE THE EAETH. 

9, Preliminary Assumptions. 

It is assumed at first that, in consequence of some internal causes, the earth is 
undergoing a deformation, but that there is no disturbance of the strata of equal density, 
and that there is no local dilatation or contraction in any part of the body. The cases 
at present excluded will be considered later. 

The result of this assumption is, that the volume of the body remains constant, and 
that the parts elevated or depressed above or below the mean surface of the ellipsoid 
have the same density as the rest of the surface. Such changes of form must, of course, 
be produced by a very small flow of the solid matter of the earth. Since the whole 
volume remains the same, this hypothesis may be conveniently called that of incom- 
pressibility ; although, if the matter of the earth flowed quite incompressibly, there 
would be some slight dislocation of the strata of equal density. 

It is immaterial for the present purpose what may be the forces which produce, and 
the nature of, this internal flow ; but it was assumed in the dynamical investigation 
that the forces were internal, and that the flow proceeded with uniform velocity. 

After deformation the body may be considered as composed of the original ellipsoid, 
together with a superposed layer of matter, which is positive in some parts and negative 
in others. The condition of constancy of volume necessitates that the total mass of this 
layer should be zero. If we take axes with the origin at the centre of the ellipsoid and 
symmetrical thereto, and let AF(0, .<p) represent the depth of the layer at the point 0, p, 
the condition of incompressibility is expressed by the integral of F(0, (p) over the surface 
of the ellipsoid being zero. Then by varying A, elevations and depressions of various 
magnitudes may be represented. 

10. Moments and Products of Inertia after Deformation. 

Before the deformation : — 

Let A, C be the principal moments of inertia of the earth ; a, b its semiaxes ; M its 
mass ; JB its mean density ; g its surface density ; and c its mean radius, so that 
oc=2a-{-b; and let the earth's centre of inertia be at the origin. 

After the deformation : — 

Let a, b, c, D, E, F be the moments and products of inertia of the above ideal shell of 
matter about the axes ; #, y, z the coordinates of the earth's centre of inertia. 

Then, since the ellipticity of the earth is small, the integrals may be taken over the 
surface of a sphere of radius c, instead of over the ellipsoid. Therefore, 

a— % 4 jjF(0, <p) sin (sin 2 sin 2 <p + cos 2 0) dd d<p 9 

Ma=hf>c 3 JJ F(0, <p) sin 2 cos <p dS dp, 

M=|7rJB^ 3 , 

and other integrals of a like nature for b, c, D, E, F, y> z. 
mdccclxxvii. 2 s 
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Since 

jjF(0, <p)$mQdd cZ<p = 0, therefore a+b+c=0. 

I£ A be the moment of inertia of the body, after deformation, about an axis parallel 

to #, through #, y, z, 

A=A+a—M(y 2 +z 2 ). 

h ~ . fh\ 2 

Now a varies as -, whilst M(y 2 +% 2 ) varies as (-] . But the greatest elevation or de- 
pression to be treated of is about two miles, whilst the mean radius c is about 
4000 miles ; hence - cannot exceed about awo? an( l accordingly the term M(y 2 +5 2 ) is 

negligible compared to a. Whence A=A-{-a. 

In like manner, the terms introduced in the other moments and products of inertia 
by the shifting of the earth's centre of inertia are negligible compared to the direct 
changes. Thus it may be supposed that the centre of inertia remains fixed at the origin, 
and that the moments and products of inertia of the earth after deformation are 
A+a, A+b, C+c, D, E, F. 

11. General Theorem ivith respect to principal Axes. 

A general theorem will now be required to determine the position of the principal 
axes after the deformation. 

Take as axes the principal axes of a body about which its moments of inertia are 
A, B, C. Let the body undergo a small deformation, which turns the principal axes 
through small angles a, j3, y about the axes of reference, and makes the new principal 
moments A, B', C. And let the moments and products about the axes of reference 
become in consequence A + a, B+b, C+c, D, E, F. Then it is required to find &, j3, y 
in terms of these last quantities. 

Let I, m, n be the direction cosines of any line through the origin, and let them 
remain unaltered by the deformation. Let I be the moment of inertia about this line 
after deformation. Let Z+W, m + lm, n + hi be the direction cosines of the line with 
respect to the new principal axes. Then, by a well-known theorem, 

ll=yryi—fin, im=ccn--yl 9 hi=fil—am. 

Now 

I = (A + a)^+(B+b)m 3 + (C+c> 2 -2Dm^-2E^-2F/m. 

But it is also equal to 

and by substituting for $1, e5m, hi, this is equal to 

A7 2 +BW + CV-2m^(a-B / )a~2fe(A'^C0i3^2/m(B f ^A f )7 
to the first order of small quantities. 
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Now this expression must be identical with the former for all values of Z, m, n ; hence 
putting Z=l, m=n = 0, A ; =A+a, and similarly B'=:B+b, 0=C-\-c. Wherefore also 

a — t ^ g, — ^ nearly, 

and 

a E F 



A-C /— B-A' 

and these are the required expressions for a, /3, y. 

If, however, B== A, y becomes infinite, and the solution is nugatory : but since, under 
this condition, all axes in the plane of xy were originally principal axes, the axes of 
reference may always be so chosen that F is zero absolutely ; and then 

D a E n 

«=— , 0=-c=A' 7=0. 



Therefore the new principal axis G is inclined to the old C at a small angle ** \ — •, 
and is displaced along the meridian, whose longitude, measured from the plane of atz 9 
is 7r+arc4an— . This is the case to be dealt with in the present problem. The posi- 
tions of the other principal axes will be of no interest. 

12. Application of preceding Theorem. 

To solve the problem numerically in any particular case, it will be necessary to find 
the integrals 

D=Ag£ 4 JJF(0, cp) sin 2 Q cos 6 sin <p dQ d<p, 

~E=h$c 4 JjF(0, <p) sin 2 6 cos 6 cos <p dftdp. 

D E 

If ^—4 and ^—4 be called d and e, then d and e stand for the above integrals, which 

depend on the distribution of surface-matter in continents and seas. 

It will be convenient to use a foot as the unit for measuring h, and seconds of arc for 
the measurement of the inclination i of the new principal axis to the old. For this 

purpose the value of the coefficient tt~t may be calculated once for all. Let its value 

when multiplied by the appropriate factors for the use of the above units be called K *. 
Now 

C-A=f( s -f)M« 2 =f(l+f s )^-f)Mc 2 . 

* I have to thank Prof. J. C. Adams for his help with respect to the numerical data, and for having dis- 
cussed several other points with me. 

2s2 
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Then if we take S—-0033439, being the mean of the values given by Colonel A. E. 

9 



Clakke, m= QQQ.a^ an( i £=20,899,917 feet*, M=|-7r|96* 3 , and — = 2, we get 



C- A=f tt^c" X '0010809 x 20,899,917 
and 

K=l-08986. 

If, in accordance with Thomson and Tait, — -=2*1, K=l - 0380, but I shall take K 

as 1*090. Then we have i' ; ==KA\/d 2 +e 2 , where K=1'090, h being measured in feet, 
and i n being the angular change in the position of the principal axis of greatest moment 
of inertia of the earth, due to a deformation given by IiF(Q, <p) all over the surface of 
the spheroid. 

The angle —7^—7— - 1 is clearly of the same order of magnitude as i, as it was assumed 

to be in Part I. 

III. FOEMS OF CONTINENTS AND SEAS WHICH PEODUCE THE MAXIMUM DEFLECTION OF 

THE POLAE AXIS. 

13. Conditions tender ivliich the Problem is treated. 

On the hypothesis of incompressibility, the effect of a deformation in deflecting the 
pole is exactly equivalent to the removal of a given quantity of matter from one part of 
the earth's surface to another. But as no continent exceeds a few thousand feet in 
average height, the removal is restricted by the condition that the hollows excavated, 
and the continents formed, shall nowhere exceed a certain depth and height. The 
areas of present continents and seas, and their heights and depths, give some idea of the 
amount of matter at disposal, as will be shown hereafter. It is interesting, therefore, to 
determine what is the greatest possible deflection of the pole which can be caused by 
the removal of given quantities of matter from one part of the earth to another, subject 
to the above condition as to height and depth. 

14. Problem in Maxima and Minima. 

This involves the following problem : — To remove a given quantity of matter from 
one part of a sphere to another, the layers excavated or piled up not being greater than 
h in thickness, so as to make \/D 2 +E 2 a maximum, the axes being so chosen as to make 
F=0. 

If D', E' be the products of inertia referred to other axes having the same origin 
and axis of z as before, it may easily be shown, from the fact that D 2 +E 2 =D' 2 +E' a , 
that D 2 +E 2 is greatest and equal to E' 2 for that distribution of matter which makes 
T) ! z=z0 and E' a maximum. 

* See Thomson and Tait, Nat, Phil. pp. 648, 651. 
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The problem is thus reduced to the following: — Rectangular axes are drawn at the 
centre of a sphere of radius c ; it is required to effect the above-described removal of 
matter, so that the product of inertia about a pair of planes through z, and inclined to 
sz at 45° on either side, shall be a maximum, subject to the above condition as to depth, 
k being small compared to c. • For convenience, I refer to the plane xy as the equator, 
to xz as prime meridian, from which longitudes -J/ are measured from x towards 3/, and 
to Q the colatitude. These must not be confused with the terrestrial equator, longitude, 
and latitude. 

A little consideration shows that the seas and continents must be of uniform depth 
&, that there must be two of each, that they must all be of the same shape, must be 
symmetrical with respect to the equator, and that the continents must be symmetrical 
with respect to the prime meridian, and the seas with respect to meridians 90° 
and 270°. 

Also the total product of inertia P, produced by this distribution, is 16 times that 
produced by the part of one continent lying in the positive octant of space ; and the 
mass of matter removed is 8 times the mass of this same portion of one continent. 

The problem is, therefore, to find the outline of the continent, so that P may be a 
maximum, subject to the condition that the mass is given. 

Take the surface-density of the sphere as unity, and let the mass removed be given 
as an elevation of a height k oyer a fraction q of the whole sphere's surface ; so that the 
mass removed from hollow to continent is Avr^kq. Then it may easily be shown that 




and 



7T 

ike* \ 2 sin 3 # sin 2-^dS 

Jo 



1 ci 

q=- ) 2>J/sinfld0, 



where \f/ is a function of d to be determined. Then writing a for 2\}/, and (Jb for cos #, 
we have to make 



j: 



{(1— ^ 2 )sina;— a cos 2 #}d^ a maximum. 



for it will be seen later that — c 2 cos 2 # is a proper form for the constant, to be introduced 
according to the principles of the Calculus of Variations. This leads at once to 

(1 — l& 2 ) cos <a=cos 2 a 
or 

sin 2 # cos 2\f/=cos 2 a, 

That is to say, the outline of the continent is the sphero-conic formed by the inter- 
section with the sphere of the cone, whose cartesian equation is 

y\l + cos 2 a) + z 2 cos 2 c6=x 2 sin 2 a. 

Eeverting to the expressions for P and £, altering the variable of integration, and the 
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limits, so as to exclude the imaginary parts of the integrals, we have as the equation to 
find & 



and 



1 r a 

*" Jo 



cos y arc cos 



/cos a \ 2 , 
Vcos x) *' 

P— 4zkc 4 I cos^\/cos 4 ^--cos 4 a^. 



These integrals are reducible to elliptic functions ; but in order not to interrupt the 
argument, I give the reduction in Appendix B. If cos 2y=cos 2 a, the result is that 

^=^ C ^?[n i (-2sinV)-F 1 ] 



cos 7 



or 



2 



and 



2=1-- {ET-F(F-E)} 

jp==. — cos y [E 1 - cos 2yF 1 ] : 



where the modulus of the complete functions E 1 , F 1 , IT 1 is tan y, and where E, F have 



a modulus _— and an amplitude 5 — y. 



cosy A " 2 

It will be observed that a is the semi-length of the continent in latitude, and y the 
semi-breadth in longitude. 

From these expressions I have constructed the following Table: — 



Semi-breadth of 


Semi-length of 


Fraction of surface 


Product of inertia 


continent 

(r). 


continent 


elevated or depressed 
to 


Wv c 


O 




/ 



•0000 


•0000 


5 


7 5 


•0054 


•0672 


10 


14 13 


•0216 


•2628 


15 


21 28 


•0486 


•5697 


20 


28 55 


•0867 


•9603 


25 


36 42 


•1362 


1-3981 


30 


45 


•1979 


1-8399 


35 


54 12 


•2732 


2*2371 


40 


65 22 






45 


90 


•5000 


2-6667 



15. Application of preceding problem to the case of the Earth. 

In the application to the case of the earth, what has been called, for brevity, the 
equator (EE in fig. 2) must be taken as a great circle, passing through a point in 
terrestrial latitude 45°. 

Figure 2 gives the stereographic projection of the forms of continents and seas 
the firm lines showing continents, and the broken ones seas, when covering various 
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fractions of the whole surface ; a and y are indicated on one of the continents. The 

other hemisphere is the same as this figure, 

when seen in a looking-glass. It will be 

observed that the limiting case is when the 

two continents fill up two quarters of the 

earth, and the two seas the other two. 

It is clear that the greatest deflection 
of the polar axis which can be produced 
by the elevation of continents of height Jc 
and having a total area 4«r6^, and the de- 
pression of similar seas, will be PK. 

A numerical Table of results will be given 
below, formed by interpolation in the above 
Table. 




IV. ON GEOLOGICAL CHANGES ON THE EAETH'S SUEEACK 

16. The points to be considered, 

It is now necessary to consider what kind and amount of superficial changes are 
brought about in the earth's shape by such geological changes as are believed to have 
taken- place. The points to be determined are : — 

i. Over what extent of the earth's surface is there evidence of consentaneous subsi- 
dence, or upheaval, during any one period. 

ii. What is the extreme vertical amount of that subsidence or upheaval. 

iii. How the sea affects the local excesses and deficiencies of matter on the earth's 
surface. 

iv. How marine and aerial erosion affect the distribution of the excess or deficiency 
of matter. 

v. The possibility of wide-spread deformations of the earth, which approximately 
carry the level surfaces with them. 

The object of this discussion is to find what areas and amounts of elevation and 
subsidence on a sealess and rainless globe are equivalent, as far as moving the principal 
axis, to those w r hich obtain on the earth. These areas and effects will be referred to as 
" effective areas and amounts of elevation or subsidence." 

It is probable that during the elevation or subsidence of any large area, the change 
proceeds at unequal rates in different parts ; probably one part falls or rises more quickly 
than another, and then the latter gains on the former. But it has been shown, in the 
dynamical part of this paper, that the axis of rotation sensibly follows the axis of figure. 
Hence it is immaterial by what course the earth changes its configuration, provided the 
changes do not proceed by large impulses, a supposition which may be certainly excluded. 
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The essential point is, to compare the final and initial distributions of matter, after and 
before a period of large geographical change. 

17. Areas of subsidence and elevation. 

When a new continent is being raised above the sea, there is no certainty as to the 
extent to which areas in the adjoining seas partake in the elevation ; even in the case 
of S. America, where the area of elevation is supposed to be abruptedly limited towards 
the west, the line of 15,000 feet depth lies a long way from the coast. 

As soon, moreover, as the land is raised above the sea, the rivers begin washing away 
its surface, and the sea eating into its coasts. The materials of the land are carried 
away, and deposited in the surrounding seas. Thus to form a continent of 1000 feet in 
height, perhaps entails an elevation of the surface of from 3000 to 4000 feet, and all 
the matter of the additional 2000 to 3000 feet is deposited in the sea. This tends to 
make the adjoining seas shallower, and to cause some increase to the area of the land. 
Therefore in a sealess globe the effect must be represented by a greater area of elevation 
and a less height. 

The bed of a deep sea is hardly at all subject to erosion, and therefore the tendency 
seems to be to make the negative features of an ocean-bottom more pronounced than the 
positive features of mountain-ranges, at least in the parts very remote from land. 

The areas, then, of existing continents may not be a due measure of the areas of 
effective elevation ; we can only say that the latter may considerably exceed the former. 
The direct evidence as to the extent of the earth's surface over which there has been a 
general movement during any one period, is also very meagre. It appears certain that 
very large portions of S. America have undergone a general upward movement within 
a recent geological period ; but there is no certainty whatever as to the limits of this 
area, nor as to whether the beds of the adjoining seas have partaken to any extent of 
this general movement. Thus the case of S. America is of scarcely any avail in deter- 
mining the point in question. The presence of deep ocean up to the Chilian coast 
seems, however, to make it probable that areas of elevation are more or less abruptly 
divided from those of rest or subsidence. 

There is only one area of large extent in which we possess fairly well-marked evidence 
of a general subsidence ; and this is the area embracing the Coral islands of the Pacific 
Ocean. The evidence is derived from the structure of the Coral islands, and is confirmed 
in certain points by the geographical distribution of plants and animals. Some natu- 
ralists are of opinion that there is evidence of the existence of a previous continent ; 
others (and amongst them my father, Mr. Charles Darwin) that there . existed there 
an archipelago of islands. In this dearth of precise information, only a rough estimate 
of area is possible. 

My father, who has especially attended to the subject of the subsidence of the Pacific 
islands, has marked for me, on the map given in his work on Coral Reefs, a large area 
which he believes to have undergone a general subsiding motion. This area runs in a 



GEOLOGICAL CHANGES ON THE EAKTH'S AXIS OF EOTATION. 295 

great band from the Low Archipelago to the Caroline Islands, and embraces the greater 
number of the islands coloured dark-blue in his map. The boundary may be defined 
as passing through :— 

N. S. 

t ' * «> t -*■ > 

Lat 3 _ j_5 _ [ 15 | 22 1 18 1 1 | 5 | 5 | 15 | 35 j 30 | 18 | 15 j 10 j 8 

Long 150] 140 | 150 | 165 | 180 | 165 j 150 | 135 | 120 | 120 | 135 | 150 | 165 | 180 | 165 

E. W. 

He also marked a smaller area, embracing New Caledonia, the S.E. corner of New 
Guinea, and the N.E. coast of Australia. 

It is noteworthy that the former large area consists of sea more than 15,000 feet 
deep, except in patches round some of the islands, where it appears to be from 10,000 
to 15,000 feet deep*. 

I marked these areas on a globe, and cut out a number of pieces of paper to fit them, 
and then weighed them. By this method I determined that the former area was "055 
of the whole surface of the globe, and the latter was "01 ; the two together were there- 
fore '065. 

It thus appears that we have some evidence of an area of between 5 and 7 per cent, 
of the globe having undergone a general motion of subsidence within a late geological 
period. But between this area and the coast of S. America there is a vast and deep 
ocean, and nothing whatever is known with respect to the movements of its bed. Hence 
it is quite possible that the area which has really sunk, in this quarter of the globe, is 
considerably larger than the one above spoken of. 

On the whole, then, perhaps from "05 to "1 of the whole surface may at various times 
have partaken of a consentaneous movement, so as to convert deep sea into land, and 
vice versa. 

Besides this kind of general movement, there have certainly been many more or less 
local rises and falls, but this small oscillation is not fitted to produce any sensible effect 
on the position of the earth's axis. 

18. Amount of Elevation, and the effects of Water. 

Humboldt has shown that the mean height of the present continents is a little less 
than 1100 feet from the sea-levelf . But this, of course, does not give the limit to the 
amount of change of level. On the other hand, there are perhaps 50,000 to 80,000 feet 
of superposed strata at most places on the earth ; but neither does this give the indi- 
cation required, because the surface must have risen and fallen many times during the 
deposition of these strata. 

But, as before pointed out, the actual upward or downward movement of land is by 

* See frontispiece-map to "Wallace's Geogr. Distrib. of Animals. 

f Sir J. Herschel seems to have doubled the height through a misconception of Humboldt's meaning. 
The mean height of the land is in English feet : Europe, 671 ; !N". America, 748 ; Asia, 1132 ; S. America, 1151. 
See a letter to ' Nature,' by Mr. J. Carrick Moore, April 18th, 1872. 

MDCCCLXXVII. 2 T 
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no means the same as its effective elevation or subsidence ; for erosion causes the effective 
to be far slower than the actual. And the actual upward or downward movement of 
an ocean-bed is different from the effective ; for the sea-water will flow off or in from 
the adjoining seas. The specific gravity of w T ater is about one third of that of surface 
rock, and the local loss or gain of matter is the actual loss or gain of surface rock, less 
the mass of the sea-water admitted or displaced. Thus the effective downward or 
upward movement of a sea-bed is about ^ of the actual ; of this a more accurate esti- 
mate will be given presently. 

It is fortunately not important to track the series of changes through their course ; 
and in order to avoid the complication of doing so, the way seems to be to estimate the 
amount of transference of matter entailed in the conversion of a deep ocean into a con- 
tinent of the present mean height. 

Suppose, then, that an ocean area of 15,000 feet in depth were gradually elevated, 
and that the final result, notwithstanding erosion, were a continent of 1100 feet in 
height. Conceive a prism, the area of whose section is unity, running vertically 
upwards from what was initially the ocean-bed. Initially this prism contained 15,000 
feet of sea-water, and finally it contains 16,100 feet of rock; so that the local gain of 
matter, on this unit of area of the earth's surface, is the difference between the masses 
of this prism, initially and finally. 

Now 1*02 is the specific gravity of sea-water, and 2*75 that of surface rock; therefore 
the same local gain of matter, in a sealess globe, would be given by an elevation of 

16,100— Iff of 15,000 = 10,436 feet. 

That is to say, 10,436 feet has been the effective elevation. 

I therefore adopt 10,000 feet as the effective elevation equivalent to the conversion 
of deep ocean into a continent ; and in the examples given hereafter, where I find the 
deflection of the pole for various forms and sizes of continent, I shall give the results of 
such an assumed conversion. 

19. Widespread Deformations of the Earth. 

It has hitherto been assumed that the elevation of land would not affect the sea-level ; 
but there can be no doubt but that elevations, such as those already spoken of, would 
do so to the extent of, say, a hundred feet. In so far, then, as this is the case, the eleva- 
tion would be masked from the eyes of geologists. But if the change of form were a 
gradual rising over a very wide area, the level surfaces would approximately follow the 
form of the rocky surface. For instance, the elliptical form of the equator carries the 
ocean level with it ; the amount of this ellipticity is such that the difference between 
the longest and shortest equatorial radii is 6378 feet *. So long, however, as these bulges 
remain equatorial they cannot affect the position of the principal axis, even should they 
vary in amount from time to time. But this kind of deformation, if not symmetrical 

* Thomson and Tait, Nat Phil. p. 648. 
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with respect to the equator, would alter the position of the principal axis, without 
leaving any trace whatever of elevation or depression for geologists to discover. 

The discrepancy which is found between the ellipticity of the earth, as deduced from 
various arcs of meridian, is, I presume, attributable to real inequalities in the earth's 
form, and not entirely to errors of observation and to the elliptical form of the equatorial 
section. It seems, moreover, quite possible that these wide-spread inequalities may have 
varied from time to time. 

Hence, even if the deposit of strata in the sea did not produce a continual shifting of 
the weights on the earth's surface, and even if geologists should ultimately come to the 
conclusion that there has never been any consentaneous elevation and depression of very 
large continents relative to the sea-level, but that the oscillations of level have always 
been local, it would by no means follow that the earth's axis has remained geographi- 
cally fixed. 

Y. NUMERICAL APPLICATION TO THE CASE OE THE EARTH. 

20. Continents and Seas of Maximum Effect. 

As far as I can learn, geologists are not of opinion that there is any more reason why 
upheavals and subsidences should take place at one part of the earth's surface than at 
another. It is accordingly of interest to suppose the elevations and depressions to take 
place in the most favourable places for shifting the axis of figure. The area over which 
a consentaneous change may take place is also a matter of opinion. 

The theorem in maxima and minima in Part III. makes it easy to construct a table 
from which that area may be selected which seems most probable to geologists. The 
following Table is formed by interpolation in the Table in sec. 14 ; the first column gives 
the fraction of the earth's surface over which an elevation is supposed to take place, a 
depression over an equal area taking place simultaneously. The second column gives 
the angular shift in the earth's axis of figure, due to 10,000 feet of effective elevation; 
as was shown in Part IV., this would convert a deep ocean into a continent. If 
10,000 feet be thought too high an estimate, the last column may be reduced in any 
desired proportion. Lastly, fig. 2 shows the forms of these continents and seas of 
maximum effect. 



Area of elevation 

or subsidence, as 

fraction of Earth's 

surface. 


Deflection of pole 

for 10,000 feet 
effective elevation. 


•001 

•005 

•01 

•05 

•1 

•15 

•2 

•5 


2% 

3 , 

22% 

1°46J' 
3° 1 f 
4° 33f 
5° 36| ; 
8° 4 J' 



N.B. The area of Africa is about -059, and of S. America about -033 of the Earth's surface. 

2t2 
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21. Examples of other forms of Continent. 

I will now apply the preceding work to a few cases where the continents and seas do 
not satisfy the condition of giving the maximum effect. 

Figures 3, 4, 5, and 6 represent the shapes of the continents as projected stereo- 
graphically. The shaded parts represent areas of elevation, the dotted parts those of 
depression ; and in the shelving continents and seas the contour lines are roughly indi- 
cated. P' shows the new position of the pole. In every case here given d=0 and F=0. 

Fig. 3. F($, <p)=sin 2$ cos 2<p, from 0—0 to sr, and <p — -~ ~ to +7, and zero over the 

rest of the globe. 



TV 

} A \\ sin 3 cos 2 cos 2<p cos <p dQ d$ =^f s/% 



#'=KAe='5480A. 

If the effective elevation or depth in the middle of continent or sea be 10,000 feet, 
PF=1° 311'. 

This is the form of continent for which $| 2 is worked out in Appendix A. 

Fig. 4. The same shape as the last, but of uniform elevation and depression of 
10,000 feet. 



7T 7T 



e=2 



4 2 j 4 sin 2 6 cos cos <p dd d<p=% V2, 

Jo J_E 



i"=Khe=V028xh. 

PP ; =2° 51 J' when 7^=10,000 ; an extreme supposition, as the area affected is a quarter 
of the whole globe. 

err* *JT™ __, 

Fig. 5. F(0, $) = !, from = to -, and from <p = — - to +-, and — \ over the rest of 
the globe. This is equivalent to F(0, <p)=f within the above limits. 

Then »"=f X 1-028 xA= -587 xh, 

and PF=1° 38', when A=10,000 feet. 

Fig. 6. F(0, <p) = sin 23 cos 2<p, from 0=0 to -, and from <p= — - to +-, and zero over 
the rest of the globe. 



2C2 



e=2 1 I sin 3 cos 2 cos 2(p cos <p d6 d<p=^, 

Jo J- 71 " 

2 



i"=KAe=-194xA; PP=32|', when A=10,000 feet. 
On the whole, then, it appears that continents, such as those with which we have to 
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deal, are competent to produce a geographical alteration in the position of the pole of 
between one and three degrees of latitude. But all these results are obtained on what 
I have called the hypothesis of incompressibility. 



Fig. 3. 



Fig. 4. 




Fig. 5. 



Fig. 6. 





VI. HYPOTHESES OF INTERNAL CHANGES OF DENSITY ACCOMPANYING ELEYATION 

AND SUBSIDENCE. 



22. A general Shrinking of the Earth. 

It may be supposed that the earth is gradually shrinking, but that it shrinks quicker 
than the mean in some regions and slower in others. This would of course lead to 
depression and elevation above and below the mean surface in those regions. A defor- 
mation of this kind may be represented as a uniform compression of the earth, super- 
posed on changes such as those considered on the hypothesis of incompressibility. If a 
be the coefficient of contraction of volume, it is clear that the values of D and E, as 

already found, must be diminished in the proportion of 1 — -«- to unity, and C— A must 

be diminished in the like proportion. Hence the deflections of the polar axis, on this 
hypothesis, are exactly the same as those already found. This seems, perhaps, the most 
probable theory, but it is well to consider others. 

The redistribution of matter caused by the erosion of continents will clearly produce 
the same effect as deformations on the theory of incompressibility. 
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23. Changes of Internal Density producing Elevation. 

In discussing the above hypothesis, I shall confine myself to the case of the upheaval 
or subsidence being of uniform height over given areas, and shall make certain other 
special assumptions. This will considerably facilitate the analysis, and will give suffi- 
cient insight into the extent to which previous results will be modified. 

I assume, then, that the elevation of the surface is produced by a swelling of the 
strata contained between distances r 1 and r 2 from the centre of the globe and imme- 
diately under the area of elevation, and that the coefficient of cubical expansion a is 
constant throughout the intumescent portion. 

This will cause a fracture of the strata of equal density, and will produce a discon- 
tinuity such as that shown in figure 7, where the dotted circle of radius r 2 indicates the 
upper boundary of the swelling strata before their intumescence. 



But the shift of the earth's axis, caused by this kind of 
deformation, will differ insensibly from what would obtain 
if there were a more or less abrupt flexure of the strata of 
equal density at the boundaries of the intumescent volume 
and of the area of elevation. 

Suppose, as before, that h is the height to which the con- 
tinent is raised above the surface ; then we require to know 
a in terms of h. 

Before intumescence, let r, 0, <p be the coordinates of any 
point within the intumescent volume; and suppose that r 
becomes r+^, whilst 9 and <p, of course, remain constant. 

The equation of continuity is easily found to be 



Fig. 7. 




du 2u 
dr* 



r 



a. 



of which the integral is %f 2 =-r 3 +/3. 

If /3 be determined, so that when r=r 1 , ^=0, 



u= 



cc r^ 



r. 



But when r=r 2 , u=li^ the elevation of the surface; therefore 



Sh 



a 



Tc 



3? 



the required expression for a in terms of li. 

* • 

Also, before intumescence, Laplace's law of internal density held good, viz. Q 8 ^-^, 
therefore afterwards the density of the stratum distant r+u from the centre is 

Q(i- a )?5ir. 
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Now the propositions given in Part II. , as to the change in the position of the earth's 
axis, remain true here also ; and the only difference is that the products of inertia D 
and E must now be expressed by different integrals. 

After intumescence the earth may be conceived to consist of \— first, itself as it was 

before ; secondly > of negative matter, of which the law of density is Q ^-^, throughout 

the space bounded by r=r x , r=c, and the cone of elevation; and, thirdly, of the 
matter which formerly lay within this space, in the configuration attained by it after 
intumescence. 

The first part clearly contributes nothing to D and E ; and the second contributes 



Q j I j r 3 sin qr sin 2 8 cos 8 j dr d8 dtp, 



integrated throughout the above space, that is from r=r x to r=c, and throughout the 
cone of elevation. 

As to the third part, the mass of any element remains unchanged, whilst its distance 
from the centre has become r+u. Hence the third part contributes 

Q 1 1 1 r(r +u) 2 sin qr sin 2 8 cos 6'- dr dd d$, 

integrated throughout the above space. 

Therefore, taking all together, and treating u as small, 



E 

Therefore 

D__E 

d e 



2Q 111 ur 2 sin qr sin 2 8 cos 81 1 ^ dr dd d<p. 
JJJ l cos <p 



3 f 

-=2Q 1 ur 2 sin qr dr, 



where d and e have the same meanings as before, in Part II. sec. 12. 

Now this last integral divides itself into two parts: first, from r=c to r=r 2 , u=.h ; 



and, secondly, from r=r 2 to r=r l9 u~^ . — g-K 

Therefore 

D E 



= 2 QJi I r 2 sin qr dr + §Q& I (r 3 — r 3 ) sin qr dr. 

Jr 2 Jr. 



u e _ j Ti 



If the value of a, be substituted, and the integrations effected, it will be found that 

D E 1 ■ 2 S , 3c 1 1 T 

.+7Z^~Z.+V~ • 7I3" 2 ~/^~\~~3~^r~ = v suppose, 



2dQAc 8 2eQhc 3 qc cos qc (qc) 2 sin qc r 2 * (^)2 /?jV sin qc 



7 \ 



r 
where S stands for the expression - sin qr + cos qr, taken between the limits c and r 2 , 

(r\ 2 r 1 2 

- ) sin qr + 2 - . — cos qr—, — w sin qr, taken between the limits 

r 2 and r v 
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Substituting in the expression i=-V—Jl — ? and using the coefficient K, we get 



i u =2KVh s /d 2 +e 2 . 

It must be noticed that this investigation is applicable as much to subsidence caused 
by internal compression as it is to elevation ; and the word intumescence is used to cover 
both phenomena. In the case of subsidence h is negative. 

Now on the hypothesis of incompressibility it was shown that ^~KA\/d 2 +e 2 . Hence, 
on the present hypotheses, the estimated deflection of the pole must be diminished in 
the proportion of 2U : 1. 

Taking £0=141° (which makes — =2, very nearly), I have calculated the values of 

2U, when- i = f M, and — =^ -^ f, i 9 0. If the earth's radius be taken as 4000 miles, 
c c 

this gives, that the superficial strata for 10 miles in thickness do not swell, but are 
merely heaved up, and that the lower surface of the intumescent volume is at the 
various distances from the earth's surface given in the first column of the following 
Table. The second column gives 2U, or the factor by which previous results would 
have to be diminished on the present hypothesis. The third column gives the so dimi- 
nished value of 1° of deflection of the pole. 



Depth below surface 

of bottom of 

intumescent volume, 

in miles. 


Factor of 

diminution of 

former results, 

(2TJ). 


A deflection of 

1 ° would be 

reduced to 

(2Uxl°). 


50 

400 

1000 

2000 

4000 


•0126 
•1011 
•2731 

•5171 
•6721 


46" 

6' 4" 

16' 23" 

31' 2" 

40' 20" 



The last row, of course, indicates that the intumescence extends quite down to the 
centre of the earth. 

This Table shows that if elevation is due to the swelling of strata at all near the sur- 
face, the alteration in the position of the polar axis would be reduced to quite an insig- 
nificant amount. The alleged deficiency of density under the Himalayas affords some 
slight evidence that it is so, at least occasionally. I believe, also, that Mr. Mallet is 
of opinion that the centre of disturbance of earthquake-shocks is not at a greater distance 
than 30 miles below the surface*. It does not, of course, follow from this evidence that 
there may not be elevations of both kinds going on, some being approximately super- 
ficial phenomena, and others probably due to unequal shrinking of the earth as a whole. 
The latter kind would be likely to produce more extensive deviations from the external 
form of equilibrium than the former. 

* Referred to at second hand by Mr. Cakrtjthers, Trans. New-Zeal. Inst. vol. viii. p. 363, 
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On the whole, then, it appears that the deflection of the polar axis cannot exceed 
that which was found in the case of incompressibility, and it may possibly be consider- 
ably less. The complete want of knowledge of the internal movements only allows us 
to state a superior limit to the change which might be produced by any one upheaval 
or subsidence. 

VII. SUMMARY AND CONCLUSION. 

24. Summary. 

For the sake of those who do not read mathematics, I will shortly recapitulate the 
chief results arrived at. 

The change in the obliquity of the ecliptic caused by any gradual deformation of the 
earth's shape of small amount is very small. Even so great a redistribution of weights 
on the earth's surface as is entailed by immense polar ice-caps during the Glacial Period, 
cannot have altered the obliquity by so much as 2^ of a second of arc ; and this is the, 
most favourable redistribution of weights for producing this effect. Thus throughout 
geological history the obliquity of the ecliptic must have remained sensibly constant. 
And, further, when the earth undergoes any such deformation, the axis of rotation 
follows, and remains sensibly coincident with the principal axis of figure. 

It thus only remains to consider the change in the geographical position of the poles 
caused by the deformation. 

The principal axes at the centre of inertia of a body are three lines mutually perpen- 
dicular, and their position is entirely determined by the shape of the body. Hence if 
a nearly spherical body be slightly deformed, the extremities of these principal axes will 
move from their original positions and describe paths on the surface of the body, which 
may be shortly described as the paths of the principal axes. In the case of the earth, 
as geologically deformed, it is only of interest to consider the path of one of these axes, 
which is, in common parlance, the earth's axis. 

If the earth be sensibly rigid, or should only readjust itself to an approximate form 
of equilibrium at long intervals (as maintained in Part I.), the geographical path of the 
axis is very nearly the same as is due merely to the geological deformation of the earth's 
shape ; but if the earth be more or less plastic, or should readjust itself frequently to an 
approximate form of equilibrium, the dynamical reactions introduced are such as more 
or less to modify the geographical path of the axis. In the case of great plasticity these 
reactions would suffice to entirely alter the character of the path. It seems probable 
that during the consolidation of the earth there was great instability in the geographical 
position of the poles. Throughout the rest of the investigation suppositions of plasticity 
are set aside, and the hypothesis of sensible rigidity is adhered to. 

Formulae for the change in the geographical position of the pole due to any small 
Reformation are found in Part II. 

On the assumption that the internal density of the earth remains unchanged by the 
mdccclxxvii. 2 u 
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deformation, the forms of continent and depression which, produce the greatest deflection 
of the poles, for the transport of a given quantity of matter from one part of the earth's 
surface to another, are then investigated. These forms are shown, projected stereo- 
graphically, in fig. 2 (p. 293). 

Part IV. gives what evidence I have been able to collect of the areas and amounts of 
deformation to which the earth may have been subjected in geological history; but as 
the discussion is not mathematical, it seems unnecessary to give an abstract thereof. 

Part V. gives numerical applications of the preceding theorems to the case of the earth, 
on the assumption that the internal density is unaltered by the deformation. From this 
it appears that the poles may have been deflected from 1° to 3° in any one geological 
period ; but the reader is referred back to that part for details. 

If upheaval and subsidence of the surface are due to a shrinking of the earth as a 
whole, but to a more rapid shrinking in some regions than others, the deflection of the 
poles is the same as that found where there is no disturbance of the strata of equal 
density. 

But if the upheaval and subsidence are due to local intumescence and contraction of 
the strata underneath the rising or falling areas, the previous numerical estimates must 
be largely reduced ; for the extent of this reduction the reader is referred to the Table 
in Section 23 (p. 302). 

It thus appears that the deflection of the poles first given is a superior limit to that 
which is possible. 

25. Conclusion. 

There remain, in conclusion, one or two miscellaneous points to be referred to. 

In a letter to Sir C. Lyell read before the Geological Society*, Sir John Herschel 
has pointed out that the isothermal strata near the surface of the earth must approxi- 
mately follow the solid surface. Therefore, when a thick stratum is deposited at the 
bottom of the ocean, the primitive bottom is gradually warmed and expands. There 
is thus a tendency for the upheaval of sea-beds, on which a large amount of matter has 
been deposited ; but this kind of upheaval certainly falls within the case of superficial 
intumescence, and could therefore affect the geographical position of the poles but little 
more than would be due merely to the weight of the deposited stratum. It must be 
noticed, moreover, that the weight of the deposited stratum would tend to compress the 
primitive sea-bed, and might counteract the expansion due to rise of temperature. 

If the earth were absolutely rigid the pole could never have wandered more than 
from 1° to 3° from its primitive position, whatever geological changes were successively 
to take place ; because the new pole could never be brought to a greater distance from 
its original position, by any fresh distribution of the matter forming the continents, than 
the maximum for this amount of matter arranged in continents of a like height. 

But it was maintained in Part I. that from time to time the earth makes a kind of 

* Proc. Geol. Soc. vol. ii. p. 549. 
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rough adjustment to a figure of equilibrium. If this adjustment is, as seems probable, 
by an earthquake, it will take place with reference to the axis of rotation at the instant 
of the earthquake. Now there exists in erosion and marine deposits a cause of terres- 
trial deformation which is certainly independent of such adjustments; and it seems 
probable that the causes of geological upheaval and subsidence are so also. We have 
therefore clearly a state of things in which the pole may wander indefinitely from its 
primitive position. On this hypothesis, as in successive periods the continents have 
risen and fallen, the pole may have worked its way, in a devious course, some 10° or 15° 
away from its geographical position at consolidation, or may have made an excursion 
of smaller amount and have returned to near its old position. May not the Glacial 
Period, then, have been only apparently a period of great cold ? If at that period the 
N. pole stood somewhere w 7 here Greenland now stands, would not the whole of Europe 
and a large part of N. America have been glaciated % And if the N. pole retreated to 
its present position, would it not leave behind it the appearance of a very cold climate 
having prevailed in those regions % 

But although such a cumulative effect is possible with respect to the geographical 
position of the pole, none such is possible with respect to the obliquity of the ecliptic. 

Now this kind of wandering of the poles would of course require extensive and 
numerous deformations, and it is hard to see how there can have been a shifting of the 
surface weights sufficient to produce it, without frequent changes in the geographical 
distribution of land and water. If, then, geologists are right in supposing that where 
the continents now stand they have always stood, would it not be almost necessary to 
give up any hypothesis which involved a very wide excursion of the poles'? 

Appendix A. (See p. 279.) 

To calculate ^ and H 2 in a supposed case of elevation and subsidence. 
Take the case of sec. 21 (fig. 3), where the elevation is given by lit sin 2d cos 2<p, from 

$=0 to sr, and from <p= — - to j, and zero over the rest of the sphere. Suppose that 

the internal motion is entirely confined to the quarter of the sphere defined by the 
above limits of 6 and <p, that radial particles are always radial, and that the motion is 
entirely meridional. 

Let #+S be the disturbed colatitude of the point #, <p. Then the equation of conti- 
nuity, which expresses that the volume of the elementary pyramid j^c 3 sin 6 dQ d<p remains 
constant, w T hen 6 becomes 0+&, is 

~tr (3- sin Q) + — sin 9 sin 29 cos 2<p = 0, 

the integral of which is 



2k 
§ sin 0+t t cos 2<p sin 3 #=a constant; 

2 u 2 
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and since S- is zero, when <p= +-, for all values of t, §= — ~t cos 2p sin 2 #, and 

-— — — ----- cos zp sm 2 9. 
Hence II 2 , twice the area conserved on the plane of xz, is 

§r 2 sin d 6?r 6?# 6?£> . — n- cos <p, 

taken from r=0 to c, Q = to 7r, <p = — - to +-. 
If the sphere be taken as homogeneous, 

H 2 = — -- § 1 1 1 r 4 sin 3 cos 2<p cos <p c?r <#S dq> 

jSTj and iZg are both clearly zero. 

The above value of H 2 is larger than what it would be in the case of the earth, if 
Laplace's law of internal density were true, because the external layers have been taken 
too heavy, and the internal too light. But taking that law of density, A=|-M<? 2 very 
nearly. 

Hence S=-i^?* 

A 5?r c 

If we let the time run on until the highest point of the continent has risen one foot, 
so that «=!__, then H i= 4 V * X 



c 20,900,000' A 5?r 20,900,000 

But reference to sec. 21 (fig. 3) shows that i"='5480/&, or in the present notation, 

&= ' 5X 64^000 nearl y- 
Therefore 

f| = _«*My*J_ nearly. 

AjS 104,500?r 2 141 J 

But generally, since the angular velocities a, j3, y of the moving axes, to which 
S?i5 H2? i? 3 refer, are very small, therefore 

7& — 77 an ir sk — jr 

to the first order of small quantities, within the limited period to which the investigation 
applies. So that in this particular case, 

1 



■—?= — -— nearly, and 5& 1 =fe 3 ==0. 

A/3 141 J ' *" ^- 3 

And, besides, this value of ^ is larger than it ought to be, because $| 2 was calculated 

on an assumed homogeneity of the earth. This, then, justifies the conclusion in the 
text on p. 279. 
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In the elevation and subsidence given by lit sin 26 sin 2<p from 0=0 to -, and from 

q> = —- to -, II x and H 2 are clearly zero, under a like supposition as to the nature of the 

— * jj 

internal motions accompanying upheavals. 



Appendix B. (See p. 292.) 

/** cos 2 a C a / — 

To reduce the integrals \ cosy arc cos — $- d% and I cos%\/cos 4 r— cos* ccd% to 

Jo cos X Jo 

elliptic functions. 



Call the former A and the latter B. 
Then integrating A by parts, 

COS 2 OL 

cos 2 ^ 
Put o?= sin %, and cos 2y=cos 2 a, then we get 



a C a • 7/ cos 2 

A=-—\ sin r a I arc cos — 5- 

Jo * \ cos ; 



* o o r^ 2sin y/ 1 -A dx 

A=2 eos2yl ( = 5—1 1 / ~ 1 — ^^==^, 

Jo \1— ^ 2 / W — 2a? 2 + sm 2 2y 



and if x=\/'l sin y sin <p, this becomes 



^/2 — ^-^{n^— 2 sin 2 y)— ■ F 1 }, where the modulus is tan y. 



Again, integrating J? by parts, 



jy—i .:•,.. 4cos 3 % sin % J % 



-J sin%. 



-1 



2 Vcos 4 ^— cos 4 a 

V 2sin y ( 1 — #2 J #2^ 

' , . — * 

V a? 4 — ■ 2.Z' 2 + sin 2 27 



T > , 7> • 1 f V2.my sm 2 2y— a? 2 — (1— a?*)a 2 _ , xT _ r x . , . , 

But i/ is also = 1 — -^======^===4-- <£# from the expression before partial mte- 

J V ^ 4 — 2^ 2 + sm 2 27 r r 

gration. Multiplying the latter expression by 2 and adding to the former, 

o -n nC ^~ sin y si^ 2 27 —-a? 2 , 
J yW 2 2y — 2# 2 + 3* 

and substituting the above value for #, 



»2 



|jg = ™?3L p+y/2 cos y(E* -- F 1 ), 

V 2 cosy v /v y 

2 \/2 
-0 = — ^— cos ypE 1 — cos2y F 1 ], the modulus being tan y. 

J? may be calculated from this form by means of the tables in Legkendke's ' Fonctions 
Elliptiques,' torn. ii. But A is not yet in a form adapted for numerical calculation. 
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The parameter — 2 sin 2 y of II 1 is negative and numerically greater than the square 
of the modulus ; therefore II 1 falls within Legendke's second class (op. cit. torn. i. p. 72). 
Now it is shown by Legendre (torn. i. p. 138) that 

^W^W*' c)-P(o)]=I+F(c)F(*, *)-E'(<?)F(J, *)-F(o)E(*, tf). 



In this case will be found to be ~— y, — ?rr-iST~ = '2 n/2 . -, and h 



cos a 



2 " A (M) cos 7 ' cos 7 

whence ^=^r-2{E 1 F~F 1 (F-E)} ? 

where the moduli of F and E are — — , and their amplitude - — y. 

cos 7 2 

From this form -4 may be calculated numerically. 



Appendix C. (Added April 1877.) 

Sir William Thomson, who was one of the referees requested by the Eoyal Society to 
report on this paper, has remarked that the subject of Part I. may also be treated in 
another manner. 

The following note contains his solution, but some slight alterations have been made 

in a few places. 

The axis of resultant moment of momentum remains invariable in space whatever 
change takes place in the distribution of the earth's mass ; or, in other words, the normal 
to the invariable plane is not altered by internal changes in the earth. 

Now suppose a change to take place so slowly that the moment of momentum round 
any axis of the motion of any part of the earth relatively to any other part may be 
neglected compared to the resultant moment of momentum of the whole*; or else suppose 
the change to take place by sudden starts, such as earthquakes. Then, on either suppo- 
sition (except during the critical times of the sudden changes, if any), the component 
angular velocities of the mass relatively to fixed axes, coinciding with the positions of 
its principal axes at any instant, may be written down at once from the ordinary formulae, 
in terms of the direction-cosines of the normal to the invariable plane with reference to 
these axes, and in terms of the moments of inertia round them, which are supposed to 
be known. 

Hence we find immediately the angular velocity and direction of the motion of that 
line of particles of the solid which at any instant coincides with the normal to the 
invariable plane at the origin. This is equal and opposite to the angular velocity with 
which we see the normal to the invariable plane travelling through the solid, if we, 
moving w 7 ith the solid, look upon the solid as fixed. Let, at any instant, x, y, z be the 
direction-cosines of the normal to the invariable plane relatively to the principal axes ; 

* This is equivalent to neglecting pf 1? f^ 2 , pf 3 of Part I. ; by which Sir W. Thomson is of opinion that 
nothing is practically lost. 



GEOLOGICAL CHANGES ON THE EAETH'S AXIS OF ROTATION. 809 

and let A, B, C be the principal moments of inertia at that instant. Let h be the 
constant moment of momentum (or twice the area conserved on the invariable plane). 

Consider axes fixed relatively to the solid in the positions of the principal axes at any 
instant, but not moving with them, if they are being shifted in virtue of changes in the 
distribution of portions of the solid. 

The component angular velocities of the rest of the universe are, relatively to these 

axes, -r-, w, -p ; and therefore, if N be the point in which the normal to the invariable 

plane at the origin cuts a sphere of unit radius, the components parallel to these axes 
of the velocity of N relatively to them are 

(h h\ (h h\ (h h_\ # 
yz y c B j, zx y^-Q), xy ^b~a 

Now, suppose that by slow continuous erosion and deposition the positions of the 
principal axes change slowly and continuously relatively to the solid. 

Let Tzr, g, a be the components round the axes (which, of course, are always mutually 
at right angles) of the angular velocity of the actual solid relatively to an ideal solid 
moving with the principal axesf . Then the component velocities relatively to this ideal 
solid of the point of the body coinciding at any instant with N are 

% %~y a -> XG—Zi*) yw—xg; 
and the components parallel to the principal axes of the velocity of N relatively to these 

(XX (a/U (AZ — — w- t 

axes are -77, ^, -j.. Hence we nave 



dx 
~dt 



dy 



dt- 1 ^ 

^ z _ 

j t =xy 









dx du dz 

These three equations give #^+^+s j£ = 0, an d therefore they are equivalent to 

two independent equations to determine two of the three unknown quantities #, y, z as 
functions of £, the three fulfilling the condition ^ 2 +2/ 2 +^ 2 =l, and it being understood 
that st, g, a are given functions of the time. 

To apply these equations to the questions proposed as to the earth's axis, let the 
normal to the invariable plane be very nearly coincident with the axis of greatest moment 

* The angular velocity of the rest of the universe relatively to the earth being opposite to the angular velo- 
city of the earth relatively to the rest of the universe, the components of the former round the axes x, y,z are 
taken as in the negative direction, i. e. from z to y, x to z, y to x. 

f «r, p, a are the same as — a, ~/3, — y of Part I. 

% These equations are the same as those given by me in Part I. p. 277. 



tj JL \J 
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of inertia C. Let O be the point where the axis C cuts the earth's surface, and let 
OX, OY be parallel to the axes of A and B. Then 2=1 ; and if the earth's radius be 
taken as unity, x and y will be the coordinates relatively to OX, OY of the point P in 
which the normal to the invariable plane cuts the surface. 

Putting therefore 2=1 in the preceding equations, we find for the determination of 
x, y that 



dx 
dt 

dy 
dt 



■yh 



B C 



e+y*> 



or 



xh(\-£)+™ 



(XiAj . 

M+ay=u, 



> 



•octr. 



(2) 



dy 

dt 



-~bx=v. 



where 



a=K 



and 



a 



B C 



<r. 



y i A. JL 



> 



• • 



• • (3) 



— a 



U= g>, V = 'G7. 



In these equations' we are to regard a, 5, u, v as given functions of the time. 
Eliminating y, we have 



d /l dx 
dt\a dt 



+bx 



d /V 

dt 



a 



-v, 



(4) 



which is a linear equation, from which x may be found by integration ; and then, by the 
first of equations (3), 

J. / (IX \ / c* \ 

y = a\ U ~~dt)'- * ,...(&) 

If B=A, the presence of <r in the equations would merely mean that the axes of x 
and y revolve with an angular velocity c; and so we lose nothing of interest with 
reference to the terrestrial problem by supposing c=0. If, then, we take A and B 
constant, equation (4) becomes, 






-\" co X- 



du 

Tt~ av > 



where 



> 



(6) 



Of 



=ab. 



To integrate this according to the method of variation of parameters, put 



and 



dx 
It' 



#=Pcosatf-|-Qsinfitf (7) 



Pa/sin^+Q&cosa^ (8) 
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so that 



d¥ , , dQi . n 

cos ^ + -77 sm.fl/i=0. 



<# 



r e find then 



P 



Q: 



1 

CO 




=J( 



* 



du 

di 

-Ti — av) cos &>£ 6?£, 



■ow ) sin^^, J 



^ 



^ 



For the case considered in Part I., where u and v are constant 



(9) 



P=- 



-$ cos cy£ + C, Q= — -5 siiifitf+C". 



and therefore by (7) 



^: 



«v 



CO 



2 



+ C cos cv£+ 6 sin wt 



(10) 



The solution expressed in equations (5), (7), (8), (9) is convenient for discontinuous 
as well as for continuously varying and constant values of u and v. 

Consider, then, the case of u=:0 and v=0, except at certain instants when u and v 

have infinite values, so that | r udt and I vdt express the components of a single abrupt 



T' 



T' 



change in the position of the instantaneous axis ; where T and T' denote any instants 
before and after the instant of the change, but so that the interval does not include 
more than one abrupt change. 

Therefore, if t be the instant of the change 



v sin cot dt= sin cot 1 vdt 

T f Jay 



v cos cot dt = cos cot 1 vdt 

Jf 



> 



T 



(11) 



Hence the part of x depending on v vanishes at the instant immediately after the 
abrupt change when t=t . Also we have by integration by parts, 



-77 sin cot dt=u sin cot— co I u cos cot dt, ] 



du 
di 



cos cot dt=u cos cot-\-co 1 u sin cot dt 



J 



r 



(12) 



And, therefore, taking the integrals between the prescribed limits, since u=0 both 
when t=T and when t=T', we have 



du . r T 

" sin cot dt=—co cos cot 1 udt, 

Jt' 



dt 



Jdu r T 

~j7 cos cot dt= cosmcot I udt 



> 



(13) 
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Using these in (9) and (7) we find, at the instant after the abrupt change, 

Ji, _— | tvLl/U » e a a » • « 9 * » » f J- TI J 

and similarly 

(J H \Aj 1/ o • • 9 * » * » » a » j I XU 

which of course might be deduced from (8) and (5). 



